Eccentricity content of binary black hole initial data 
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Using a post-Newtonian diagnostic tool developed by Mora and Will, we examine numerically 
generated quasiequilibrium initial data sets that have been used in recently successful numerical 
evolutions of binary black holes through plunge, merger and ringdown. We show that a small but 
significant orbital eccentricity is required to match post-Newtonian and quasiequilibrium calcula- 
tions. If this proves to be a real eccentricity, it could affect the fine details of the subsequent 
numerical evolutions and the predicted gravitational waveforms. 
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I. INTRODUCTION 

Recent breakthroughs in numerical relativity have made it possible to evolve Einstein's equations for binary black 
holes (BBH) stably for several orbits, including the plunge, merger and ringdown phases, and to generate intriguingly 
robust gravitational waveforms [ESSSQ- The starting point of these evolutions is a set of initial data, obtained from 
the initial- value equations of general relativity, intended to represent two black holes in circular orbital motion; this is 
the expected end product of long-term binary evolution under the circularizing and damping effects of gravitational 
radiation emission. 

In earlier work 0, Q we developed an approach, based on the post-Newtonian approximation, designed to study 
and elucidate the physical content of these initial data sets, and showed that, in order to match post-Newtonian 
theory with some data sets ||, a small but significant orbital eccentricity was required. In this paper we apply 
this post-Newtonian diagnostic to the initial data used in recent BBH evolutions, and find that they also require an 
orbital eccentricit y. I n particular we examine the corotating and non-spinning initial data computed by Cook and 
collaborators |lCt 111) and the non-spinning "puncture" initial data of Tichy and Briigmann and of the "Lazarus" 
group fl^ . Il3j . If this residual and unintended non-circularity is real, it may affect the detailed structure of the 
numerically generated gravitational waveforms. 

The plan of the paper is as follows. In Sec. [H] we summarize the Post-Newtonian diagnostic equations for BBH 
derived in Q. In Sec. IIIII we apply the diagnostic to the BBH quasiequilibrium configurations of 0, 0, 0, Il3| . 
Sec. IIVI presents conclusions. 

II. POST-NEWTONIAN DIAGNOSTIC FOR BINARY BLACK HOLES 

Consider a binary system of black holes of irreducible masses mi and m,2, and rotational angular velocities u\ 
and a->2, with mj rr = mi + m-i and r\ — mxmilrri^^ defining the total irreducible mass and reduced mass parameter, 
respectively (0 < r/ < 1/4). Following [|| Q we define an eccentricity e and a quantity ( = m- m jp related to the 
semi-latus rectum p of the orbit, according to: 



/ \ 4/3 
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where Q p is the value of the orbital angular frequency f2 where it passes through a local maximum (pcricenter), and 
S7 a is the value of f2 where it passes through the next local minimum (apocenter). These quantities reduce exactly 
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to their Newtonian counterparts in the small orbital frequency (Newtonian) limit, and are gauge invariant through 
second-post Newtonian order, among other advantages @>0- 

We want to compare with quasiequilibrium configurations of equal-mass BH-BH binaries, so we set TOi = mi and 
77 = 1/4. For corotating binaries we also set u)% = L02 = u) = £1, while for non-spinning binaries we have u>i = L02 = 0. 
We exploit the fact that there exist exact formulae for the energy and spin of isolated Kerr black holes in terms of the 



irreducible mass, M = M irr /[1 - 4(M irr u;) 2 ] 1 / 2 , S = AM? r u/[\ 
momentum of the system are then given by 
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The orbital ("Orb") contributions are expressed in the ADM (Arnowitt-Deser-Misner) gauge and are valid to third 
post-Newtonian (3PN) order. In Eqs. (|3a|) and 13b|l . we have expanded the Kerr formulae for M and S in powers of 
TOi rr a>, assumed to be small compared to unity, keeping as many terms as needed to reach a precision comparable to 
our 3PN orbital formulae, and have subtracted TOi rr in order to obtain the binding energy. The "N,Corr" terms come 
from converting the individual total masses that appear in the Newtonian orbital energy to irreducible masses and 
their corrections due to spin, and the 1 '&pin" terms are spin-orbit effects. For black hole binaries, tidal and spin-spin 
effects can be shown to be negligible 0- To obtain Eb and J at a turning point as functions of f2, we substitute 
C = (mi rr f2 a ) 2 / 3 /(l — e) 4 / 3 or ( = (TOj rr f2 p ) 2 / 3 /(l + e) 4 / 3 for apocenter or pericenter, respectively. When E^, J, uj 
and are suitably scaled by m- m , there remains only one free parameter, the eccentricity of the orbit, This approach 
was used in to compare with the numerical quasiequilibrium solutions of Grandclement et al. , and it was found 
that a substantially better fit to the numerical data was obtained for non-zero values of e, of the order of 0.03, with 
the system at apocenter, than for e = 0. We now apply this diagnostic to other data sets that have recently played 
an important role in BBH evolutions. 



III. DIAGNOSIS OF BBH INITIAL DATA SETS 



Cook and collaborators have developed initial data sets for quasiequilibrium BBH, allowing for both corotation 
and zero spin, in a series of papers 0, 0> Q, 0, 0, 0, 0|. Using the thin-sandwich approach, combined with 
"excision" boundary conditions for the black holes adapted for treating spin, they considered systems possessing a 
"helical Killing vector" , d/dt + fid/ dip, meant to represent a circular orbit, one that is stationary in a frame rotating 
with angular velocity fl. Additionally, they impose the condition that the Komar mass, a mass defined for stationary 
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FIG. 1: Binding energy (left) and angular momentum (right) for the corotating and non-spinning BH-BH initial data from 
fTl| . Circles (squares) are non-spinning (corotating) data from [Tl|. the solid lines are circular PN diagnostics with e = 0, the 
dashed lines are eccentric PN diagnostics with e = 0.025. 




FIG. 2: Eccentricity required to match data points for non-spinning (left) and corotating data (right). Black lines are solutions 
for the energy, red lines are solutions for the angular momentum. In the non-spinning case, continuous lines refer to the Caudill 
et al. data fll| . dashed lines to the earlier Cook-Pfeiffer data fTo||. 



systems, equal the ADM mass, an invariant mass measured at spatial infinity. It is believed that this condition helps 
ensure that the orbit is truly circular. We apply our diagnostic to two data sets, taken from Refs. P5 and [Tl| . 
respectively. For non-spinning BH, the second data set used a more accurate prescription for setting the BH spins to 
zero; in the earlier data, the black holes were not truly nonrotating. We take the data from Tables IV (corotating) and 
V (non-spinning) of [Tol ] and of and plot Eb/m,i TI and J/m? r vs. mi rr £l. Figure ^ shows the comparison between 
the data of Caudill et al. an d our diagnostic, plotted for e = and for e = 0.025 (with our definitions, positive 
values of e correspond to the system being at apocenter). Figure |U shows the eccentricity required to match each 
data point from both pH | and In the improved data set of Caudill et al. for the non-spinning case, the apparent 
eccentricity in the fits to J is reduced, and the functional behavior of e with mi rr J7 is now the same (monotonically 
increasing) in both the non-spinning and corotating cases. For the corotating case, there is essentially no difference in 
the fits between the two data sets. Furthermore, as in earlier comparisons |6(|, there is a systematic difference between 
the eccentricity required to match the binding energy and that required to match the angular momentum. 

Another approach to initial quasiequilibrium data for BBH evolutions is the "puncture" method, in which the 
conformal factor of the conformally flat spatial slices is written in terms of a Newtonian-like potential m^/lx — x^| (a 
"puncture" ) for each body. This approach can also be made to incorporate the helical Killing vector and Komar-ADM 
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mass equality assumptions [lflj . Tichy and Briigmann [T^ and Baker et aL [lflj (the "Lazarus" project) have used 
this approach to generate initial data for non-spinning binary black holes in quasi-circular orbits. In Figure[31we show 
the eccentricity required to match these data sets; again the eccentricity is small but significant, and again different 
between Ef, and J. We also show explicitly the eccentricity required to match the Tichy-Brugmann data using a 2PN 
diagnostic, obtained by truncating the 3PN terms in the orbital expressions; the differences are comparable to the 
differences between the two numerical data sets, and are not large enough to account for the difference between the 
-Eh and J curves. Table [I] lists the coefficients of a cubic fitting function to the eccentricity required to match Ej, and 
J for different initial data sets. 

TABLE I: Fit of the eccentricity for the corotating (corot) and non-spinning (nospin) data by Caudill et al, for the Tichy- 
Briigmann puncture data (TB) and for the Lazarus data (Lazarus). The integer TV is the number of data points used for the fit. 
We carry out a least-squares fit by a cubic polynomial e = X^fc=o e k{rni TT 0) k . Ae max = max[(e — e num )/e llum ] is the maximum 
percentage error of the fit with respect to the numerical data e Ilum . 
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IV. CONCLUSIONS 



We have shown, using a post-Newtonian diagnostic tool, that initial data sets for binary black hole mergers may 
actually represent slightly eccentric orbits. Several remarks are called for. First, there is evidence from the dynamical 
evolutions of some of these initial data sets that the orbits are slightly eccentric. For example, in the recent evolution 
of several BBH orbits through merger and ringdown by Baker et al. [3( using the Lazarus initial data, an oscillatory 
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behavior of the separation of the black holes can be seen in their Figure 9. Similar oscillations were seen in the binary 
neutron star evolutions of |20l |21| , although the results there were very sensitive to grid resolution and size of the 
computational domain. 

On the other hand, we continue to be puzzled by the difference in values of e inferred from fits to Ef, and J. This 
difference was also seen in fits of the diagnostic to data from the Meudon group |(|, and could be cited as a defect 
of the PN approximation. However, this difference occurs systematically even at the smallest values of TOi rr f2, where 
relativistic corrections are quite small. 

We want to emphasize that the eccentricity we are discussing here is not related to the mismatch between an initial 
quasicircular orbit (with r = r = by construction) and the reality of a pre-existing inspiral (with r m — 16(m/r) 3 /5), 
since the initial data sets know nothing about radiation reaction. That eccentricity, which would be induced on an 
evolution from a perfectly circular initial orbit, has been discussed in detail by Miller J^]. Depending on the starting 
point of the evolution, the induced eccentricity from this effect could be as large as 0.03. 

Irrespective of the origin of the eccentricity, Miller pointed out that the result could be a substantial decrease 
in detection signal-to-noise when a numerically generated, eccentric waveform template is matche d ag ainst a "true" 
waveform generated by a quasicircular inspiral of a real BBH (see, for example Figures 7 and 9 of [22j). 

If eccentricity is an issue and cannot be removed or reduced by tuning the initial data sets, one could ask whether 
it could be damped away naturally by numerically evolving several orbits leading up to the onset of plunge, around 
TOi rr ri ~ 0.1. Using Eqs. (2.34) of Q, which give the evolution of our orbit elements e and £ under radiation reaction, 
it is straightforward to show, at 2.5PN order and in the small eccentricity limit, that the number of orbits N required 
to reduce the eccentricity by a factor X = e//ej by the time the orbit reaches a final angular velocity f2/ is given by 
N = (X -30 / 19 - l)/647T77(mi rr fJ/) 5 / 3 . For r] = 1/4 and m^rO/ = 0.1, this gives 34 orbits for a reduction by 1/10, 
and 11 orbits for a reduction by 1/5. Suppressing eccentricity this way is likely to be a challenge without additional 
breakthroughs in numerical relativity. 
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